Keldysh study of point-contact tunneling between superconductors 
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We revisit the problem of point-contact tunnel junctions involving one-dimensional superconduc- 
tors and present a simple scheme for computing the full current-voltage characteristics within the 
framework of the non-equilibrium Keldysh Green function formalism. We address the effects of dif- 
ferent pairing symmetries combined with magnetic fields and finite temperatures at arbitrary bias 
voltages. We discuss extensively the importance of these results for present-day experiments. In 
particular, we propose ways of measuring the effects found when the two sides of the junction have 
dissimilar superconducting gaps and when the symmetry of the superconducting states is not the 
one of spin-singlet pairing. This last point is of relevance for the study of the superconducting state 
of certain organic materials like the Bechgaard salts and, to some extent, for ruthenium compounds. 
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I. INTRODUCTION 

The theory of superconductivity by Bardeen, Cooper 
and Schrieffer (BCS) is one of the most important 
achievements of condensed matter theory. Some of the 
most striking consequences of this theory concern the 
tunneling to and from a superconductor. Indeed, the his- 
tory of tunneling experiments and applications is strongly 
linked to that of superconductivity. Not only some of 
the most crucial experimental verifications of the BCS 
theory came from tunneling experiments, A but also some 
of the most important practical applications of super- 
conductivity involve Josephson tunneling junctions. To 
describe the manifold of experimental and practical situ- 
ations, two limiting cases are usually considered: planar 
interfaces and point contacts. As of latter, point-contact 
tunneling per se acquired renewed relevance with the de- 
velopment of scanning tunneling microscopy (STM)f2*2*A 
that is today at the forefront of the experimental tech- 
niques used to study unconventional superconductors. 
For STM the tip can be modeled using some idealized ge- 
ometry. For example the cases of sphericalj££ conical^ 
and pyramidal 9 tip geometries were considered in the lit- 
erature (the last two were used to model very close STM 
contacts). The point-contact approximation is therefore 
the simplest one for the kind of tunneling processes that 
take place on STM experiments. Other ways to real- 
ize point contacts include the use of break junctions and 
pressed crossed wires. 

The simplest theoretical models used to interpret 
experiments involving superconducting tunneling are 
typically based on a simple scattering picture and 
go generally under the name of semiconducting band 
uitjilrlsM^LLLiJA A more systematic approach is that 
based on the tunneling Hamiltonianii2ii£ii£i A large se- 
ries of recent expeiimentai^^^^^^ on atomic-size 
contacts showed impressive agreement with the theory; 
achieving, some of them, detailed microscopic description 
of the contacts. The current transport in these systems 
can be described as taking place through a small num- 
ber of independent conduction channels, each of them 
well described by a point-contact model. In some exper- 
iments even the observation of single-channel transport 
was possible. 



Tunneling can thus be used as a very efficient probe 
of the properties of the leads. In particular one can 
expect to use it to determine the symmetry of the su- 
perconducting order parameter in the leads. However, 
the previous theoretical analysis of point-contact tun- 
neling, although efficient in simple cases, are too cum- 
bersome to be easily generalized to more complex cases 
such as unconventional order parameters at finite tem- 
peratures and finite magnetic fields. Simplified semi- 
classical methods exist ^^Sk but they suffer from their 
own limitations, for instance when one is dealing with 
anisotropic superconductors^ Thus a general and simple 
microscopic theory of point-contact tunneling was clearly 
lacking, and is necessary in order to take into account 
some of the complications of unusual superconductivity. 
Providing such a theory is the purpose of this work. We 
use a Keldysh formalism, to be able to compute the full 
current-voltage characteristics and gain access to the ef- 
fects of external magnetic fields, potential scattering bar- 
riers and finite temperatures on the transport properties 
of different junctions at arbitrary finite voltages. Con- 
trarily to previous implementations of this technique, us- 
ing the solution of difference equations^ we here obtain 
and diagonalize the full tunneling action for the point 
contact tunneling junctions involving normal-metal and 
superconductor leads. This allows one to easily incor- 
porate complications such as triplet pairing in the leads, 
finite temperature and finite magnetic field. We explore 
in particular the physical properties of tunneling systems 
with leads with triplet paring parameters. 

Indeed, although the possibility of having triplet pair- 
ing was investigate d^'^ ! soon after the BCS theory, and 
such an unconventional scenario was found about the 
same time in the p-wave spin-triplet superfluid state 
of 3 n e ^^LSL the quest to identify a p-wave charged 
superfluid proved much more challenging. A class 
of candidates for triplet pairing, though the evidence 
is as yet not completely conclusive, are the organic 
superconductors^^ and the ruthenates^*^ There are 
also proposals of spin-triplet pairing phases for some 
heavy fermion superconductors like UPt3, but the is- 
sue remains more open in those cases^i The organic 
compounds are the most interesting for us due to the 
quasi-one-dimensional nature of their normal phases, and 



2 



also because there is currently considerable debate on 
the symmetry of the superconducting phasem^^^ For 
the ruthenates, on the other hand, the triplet pairing 
seems already backed up by a considerable amount of 
experimental evidence 42ai To sort out this question of 
the symmetry of the order parameter, tunneling can thus 
be an invaluable tool. Recently, STM tunneling exper- 
iments were used to study the symmetry of the super- 
conducting phase of Sr 2 Ru04 and other compounds. No 
such attempts were as yet made in the case of the quasi- 
one-dimensional organic salts, but efforts in this direc- 
tion are on their way. Also recently, preliminary experi- 
ments involving junctions between two (different) Bech- 
gaard salts were performed, and they showed a number of 
puzzling features including a zero-bias conductance peak 
("anomaly") and zero excess current^ In that context, 
a precise theory of the particularities of point-contact 
tunneling involving spin-triplet superconductors, done in 
the microscopic framework of the tunneling Hamiltonian 
models, was called for. Given the nature of these systems, 
it is difficult to perform phase sensitive experiments such 
as the ones that were, for the cuprates, smoking guns 
to fix the symmetry of the order parameter. We show, 
however, that the tunneling spectrum has characteristic 
features, such as the magnetic field dependence, that can 
be used to unambiguously determine the order parame- 
ter symmetry in these systems. A short account of part 
of the results of this paper was published previously^ 

The rest of the paper is organized as follows. In 
Sec. [H] we present the model that we use to describe 
the point-contact junction geometry between either nor- 
mal or superconducting leads. In Sec. IIIII we work out a 
way of finding the tunneling characteristics using a non- 
equilibrium (Keldysh) formalism. This allows us to ob- 
tain the current- voltage characteristics for arbitrary volt- 
age, temperature or magnetic field, for junctions with ei- 
ther normal, or singlet or triplet superconducting leads. 
All technical details have been confined to these two sec- 
tions, while the two remaining sections deal with the 
physical consequences of our findings. Readers only in- 
terested in those can thus safely jump to Sec. IIVI where 
the physics of such junctions is discussed in detail. Those 
results are put in context within different experimental 
possibilities in Sec. In particular we discuss there 
the possibility of using tunneling experiments to probe 
the nature of the superconducting pairing in the organic 
superconductors. In Sec. IVII we close the paper with a 
general discussion of the implications of our results. 



II. MODEL OF THE POINT-CONTACT 
JUNCTION 

Using a non-equilibrium Keldysh formalism we calcu- 
late the full current-voltage characteristics of different 
types of tunnel junctions where each side of the junc- 
tion can be either a normal metal (N), a singlet (S) or 
a triplet (T) superconductor. We start from a tunneling 



Hamiltonian formulation, 

H = H 1 +H 2 + H tun (1) 

H tun = *«' vi(o)Vv CT (o) (2) 

t,l',<T 

The first two terms describe the two leads of the junction 
(superconducting or otherwise) and the third one models 
the tunneling processes in which an electron with spin a 
hops from lead £' into lead £. The tunneling matrix is 

The diagonal terms, V n , are local contact potential terms 
included for the sake of generality^ and the off-diagonal 
ones are the tunneling matrix elements taken to be con- 
stant consistently with the assumption of a point contact. 
Since the number of particles in each lead is a conserved 
quantity in the absence of tunneling, we can define the 
current as proportional to the rate of change in the rela- 
tive particle number and writei^ 

I=Udt (N 2 - Ni)) = £ ([JT tun , Ni - N 2 }) . (4) 
2 it 

Notice that the diagonal part of the tunneling matrix 
conserves particle numbers and will not contribute to the 
current. 

To model the superconducting leads in calculations 
intended to capture the main features of point-contact 
transport on conventional superconductors, very simple 
models suffice to achieve even quantitative agreement 
with the experiment. Contrary to the case in some pla- 
nar junction experiments, dimensionality plays little or 
no role in the tunneling. Therefore one can use one- 
dimensional leads to carry out all the standard calcula- 
tions. The situation becomes more complex in the case 
of unconventional superconductors, mainly because the 
anisotropic nature of the pair wave-function has to be 
taken into account when modeling the leads. The most 
conspicuous case is that of the cuprate compounds, for 
which the putative <i-wave paring cannot be modeled 
within a single-band one-dimensional lead. On the other 
hand, the organic superconductors that we are interested 
in are supposed to have p-wave symmetry. Since both 
s-wave and p-wave symmetries can be modeled in single- 
band one-dimensional chains, we can conveniently set up 
a formalism that encompasses the two cases, as well as 
the normal state. In the following, we will consider a one- 
dimensional band with two Fermi points and expand the 
fermion fields around them in the conventional way4^ 

^ a (x) « e~ ik ^ La (x) + e*'*^ Or) (5) 

thus defining left and right moving fields (lead indexes 
were omitted here). Using these fields and in the spirit 
of the BCS theory, we introduce the following four gap 
functions: 

A a (x) = A a (a V ia (a:) a a a0 i> Rj3 (x)) (6) 
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where Greek indexes are summed over, a = 0, . . . , 3 and 
cP a p is the identity matrix while the other three are the 
usual Pauli matrices. We use the notation a — —a with 
a £ (|,T) = The constants A a would de- 

pend on the details of the microscopic pairing mechanism 
about which we make no assumptions. 37 With this def- 
inition Ao (x) is the spin-singlet order parameter, as in 
conventional superconductors, and the other three func- 
tions form a vector of spin-triplet order parameters^ 
A (x) = A (x) d (x) . We use the approximation of drop- 
ping the spatial dependence in the order parameter and, 
directly in Fourier space, we write the Hamiltonian for 
any of the two leads as 



where K = H — fiN with /x the chemical potential of 
that lead. All the indexes are summed over, in particular 
c G (L, R) = (—1, +1) sums over the two possible chirali- 
ties and ^ c fe CT = cv-pk — fi—ah are the corresponding linear 
dispersions, shifted by the inclusion of chemical potential 
and magnetic field along the i-axis (for convenience we 
will take i>f = 1). This is the natural extension to the 
triplet case of the usual pairing-approximation Hamilto- 
nian found in BCS theory, we remark that the fact it does 
not conserve particle number is an artifact of the anoma- 
lous mean field approximation behind its derivation and 
has no bearing in the operator definition of the current. 



III. LOCAL ACTION APPROACH 



Here we arranged the different components of the order 
parameter using the following matrix notation: 



_ i Aii A n\ = a *a _ ( Ao + A 3 Ax - iA 2 
'All Afi J ~ a -U 1 + zA 2 A -A 3 



We note that another convention, the one introduced in 
the work of Balian and Werthamer^ is related to ours 
via Abw = A - (ia v ); the difference is rooted in a different 
definition of the spinor basis. 

In the case of zero magnetic field, K 2 is block diagonal 
and one arrives to a closed solution for the quasiparticlc 
excitation spectrum^ In the presence of magnetic field 
the calculations for the case of a general order parameter 
are more involved. We adopt the convention of taking 
the quantization axis ( z ) along the magnetic field direc- 
tion and consider the cases of triplet order parameters 
parallel or perpendicular to the field. In both of these 
cases the Hamiltonian can be diagonalized via a canon- 
ical rotation (i.e. a Bogoliubov-Valatin transformation) 
that proceeds in a completely identical way to that of the 
conventional s-wave case. Following the analogy further, 
the local Green functions for the leads can be written 
down immediately. For the case of a parallel order pa- 
rameter (i.e. Ai = A2 = 0) the non-zero matrix elements 
of the advanced and retarded Green functions are: 



Within the extended-BCS framework, the Hamiltonian 
remains a quadratic form including 'anomalous' terms. 
To be able to write it down as a canonical quadratic form 
we introduce the following spinor notation: 



kner 



(w) 



^Rka M 



(7) 



(that we present directly in Fourier space) . Here k is the 
reduced momentum (after linearization was carried out) 
and 



w = iv - 



(8) 



r.a r.a 

9u =533 



r.a r r,a-\ 
9l3 = [931 \ 



r.a r.a 

922 — 9aa 



r.a r r.ai 
524 = [542 \ = 



(vj + h ± irf) 



"' '\A n \ 2 -(w + h±ir,) 2 



'11 



w ^\A u \ 2 -(w + h±ir ) f 
2 — (zu — h ± it]) 



"' l \A n \ 2 -(vj-h±ir{f 



Ml 



"' '\A n \ 2 -(w-h±i V f 



(10) 

(11) 
(12) 

(13) 



is the shifted frequency corresponding to a time evolution 
given by K (cf. with the discussion of tunneling given in 
Ref. |4y|); here the bars have the meaning of minus signs. 
Since we make explicit distinction between chiralities, all 
the components of the spinor are independent. Using this 
basis the Hamiltonian can be written in matrix form, 



kner 



(w) 



ar 
.0 

(XT. 



ffcmr M (9) 



with the upper (lower) sign corresponding to the retarded 
(advanced) ones. Here w = 4up is an energy scale re- 
lated to the Fermi velocity (or equivalently to the normal 
density of states at the Fermi level) and r\ is a positive 
infinitesimal that regularizes the Green functions (some- 
times kept finite to model the inelastic relaxation pro- 
cesses inside the leads). Analogously for the case of a 
perpendicular order parameter (i.e. Aq = A3 = 0), the 
non-zero matrix elements of the advanced and retarded 
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Green functions are this time: 

2 — {vj ± irf) 



r,a r,a 



r.a r r.ai 
#14 = L941 \ 



r,a r,a 

.922 — #33 



r.a r r.ai 
923 = 1.932 J 



" , \A n f-(«,±i V j l 



TT 



(14) 

(15) 
(16) 

(17) 



The non-equilibrium formalism that we seek to imple- 
ment in order to access the full I-V characteristics for 
arbitrary finite voltages, requires the introduction of 
one more linearly independent Green function. From 
the expressions for the retarded and advanced func- 
tions and using the assumption of thermal equilibrium 
of the leads, we can construct immediately the so-called 
Keldysh component^ of the local lead Green function: 
( 5 r.-4)tanh(w/2T). 
Before proceeding, we stop to comment on the case of 
two normal leads. The corresponding Green functions are 
obtained by using any of the two sets above and taking 
the limit A a — > Va. In this case it is a simple exercise 
to derive from Eq. (0J the well-known expression for the 
conductance of an N-N junction: 



Gnn = -^a with a 
7m 



At 2 



2 • 



(18) 



where we reintroduced Plank's constant and measured t 
in units of w. This expression was first derived by Lan- 
dauer and later extended and generalized in the works 
of Biittiker, Imry and others, ~ The constant a is called 
the channel transparency and takes values in the interval 
[0, 1]. Now we return to the case when at least one of the 
two leads is superconducting. 

Given the local Green functions and the tunneling 
Hamiltonian, the simplest way to proceed in order to 
compute the characteristics of a junction is to use linear 
response and perturbation theory^£*2& A more rigorous 
approach should make use of non-equilibrium Green func- 
tions and treat the tunneling term to all orders. This, to 
be able to calculate the full I-V line and give a quantita- 
tive account of its sub-gap structure even in the ballistic 
limit (i.e. for a — > 1). One past implementation of this 
program made a clever use of the non-equilibrium Dyson 
equations and reduced the problem to the solution of a 
set of linear recursion relations. 16 Here we do not want 
to restrict ourselves to the s-wave case and to zero tem- 
perature and fields, we shall take then a different route. 
We treat the local action and Green functions directly as 
matrices in order to gain the convenience of a simpler im- 
plementation of multiband multicomponent spinors and 
deal with them numerically. 



We notice that the lead Green functions can be in- 
verted in close analytical form to obtain the correspond- 
ing local lead actions. This procedure implies the as- 
sumption of fast relaxation rates, 1 which is consistent 
with the point-contact geometry of the junction. Namely, 
working in a Keldysh-extended Nambu-Eliashberg spinor 
basis of symmetric and antisymmetric combinations of 
forward and backward time paths (see Ref. |4jJ), the lo- 
cal action for a single lead can be written as 



57 



2^ 



(19) 



where n labels the different components of the four- 
spinors as introduced in Eq. 10 and k is the index for 
the two (symmetric and antisymmetric) Keldysh compo- 
nents. The matrix representation of the spinorial action 
density is given by 



k9 9 



(- w 



] _1 9 k [sT 1 
[sT 1 



o 



(20) 



where g r ' a ' k are the matrices in the four-spinor basis 
whose nonzero components were given above (for the two 
orientations of the order parameter that we will consider, 
the inverses of g r,a are easy to calculate in closed form) . 
By combining these actions and the spinorial matrix rep- 
resentation of the tunneling Hamiltonian (7Ytun) writ- 
ten in a two-lead Keldysh-extended Nambu-Eliashberg 
spinor basis, one can ensemble the full non-equilibrium 
action matrix density for the junction 



A = [At=i © At=a] - «tv 



(21) 



While carrying out this construction, special attention 
must be paid to the fact that the shifted frequencies [see 
Eq. (JHJ] will have different reference levels when there 
is a relative bias applied to the leads. Positively- and 
negatively-shifted-frequencies in each lead are related by 
the coherent pairing processes in the superconductors; 
this is reflected on the choice of frequency pairs in the 
spinor basis. When two superconductors with different 
chemical potentials are put into contact, the tunneling 
Hamiltonian connects real (i.e. unshifted) frequencies. 
Thus, at finite voltages, pairing and tunneling together 
create an infinite set of related frequencies that is at the 
heart of the multiparticlc tunneling processes mediated 
by the so-called Andreev reflections; this is illustrated in 

pig. in 

To each value in the frequency window defined by the 
chemical potentials in the two leads, one such set of 'en- 
tangled' frequencies can be assigned. These sets are in- 
dependent and the action is block diagonal between dif- 
ferent ones. Discretizing the frequencies in this window 
automatically defines a discretization of the 'whole' fre- 
quency space. We proceed in this way and deal with one 
such set of frequencies at a time. Since these sets are 
infinite, we truncate their hierarchies at some distance 
from the central frequency window. This is equivalent 
to introducing a soft limit in the number of allowed An- 
dreev reflections: up to some fixed number (-/Va) they 
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FIG. 1: Depiction of a set of frequencies involved in a multi- 
ple coherent tunneling process (the vertical axis corresponds 
to frequencies). The horizontal lines correspond to frequency- 
conserving tunneling processes that take place with amplitude 
\t\, while the vertical lines correspond to electron (hole) pair 
creation or destruction processes that occur with amplitude 
|A| (the superconducting gap is taken in this figure to have 
the same magnitude in the two sides of the junction). The 
chemical potentials, thermal distributions and superconduct- 
ing quasiparticle densities of states in the leads are schemati- 
cally indicated. The lowest order nonzero contribution to the 
tunneling is depicted. The dashed lines and arrows indicate 
that the chain of interrelated higher order processes continues 
ad infinitum. 



are fully taken into account and then they are gradu- 
ally suppressed until twice that number is reached. This 
is a natural and consistent cut-off scheme at any finite 
voltage, since the presence of a growing frequency de- 
nominator makes the higher contributions less and less 
important regardless of the value of a. It is also clear 
what the limitations of the approach are: as the differ- 
ence in chemical potentials decreases, the denominators 
grow more and more slowly and a larger number of An- 
dreev reflections is required in order to achieve the same 
accuracy. 

The implementation of the described scheme imports 
one more complication. In the spinor basis we adopted, 
only chirality conserving tunneling processes can be writ- 
ten in matrix form. To overcome this problem we in- 
troduce a second mirrored spinor basis with the chiral- 
ities inverted (as in Eq. (JJ but interchanging R <-> L). 
Using two copies of the spinor space the full tunneling 
Hamiltonian (that is, including chirality non-conserving 
processes) can be written as a matrix and the whole fre- 
quency space for both chiral species is considered (let us 
stress that no Hilbert space doubling takes place). Invert- 
ing the action matrix density thus constructed [Eq. 1|21|) ]. 
using standard numerical methods, we can obtain fre- 
quency densities for the different current harmonics (con- 
structed out of the Keldysh components of lead-mixing 
Green functions). Here we will concentrate on the dc 



component. Finally, the current is computed integrating 
its density over the full frequency axis, 

The practical implementation of this numerical scheme 
is straightforward and allows one to consider the (com- 
bined) effects of finite temperature, applied magnetic 
fields, contact potentials in the junction, spin-flip tun- 
neling or spin-flip scattering processes^ in the leads. It 
is also possible to compute the ac response. Here our 
primary interest is in comparing singlet and triplet su- 
perconductor junctions and how they respond differently 
in the presence of an external field or with tempera- 
ture; other additional complications will be discussed 
elsewhere. Even though our numerical scheme is not 
well suited for studying the limit V — > 0, in particular 
the combination a ~ 1 with V ~ is the computation- 
ally most expensive one, that limit can nevertheless be 
approached analytically. On the other hand, all other 
regimes can be solved with modest computational effort 
and the algorithm is quite easily parallelizable. Even 
more, we shall show that the finite bias features are the 
ones that might provide useful signatures to further es- 
tablish the spin-triple scenario. 

IV. TUNNELING CHARACTERISTICS 

To discuss and compare the I-V characteristics for dif- 
ferent types of junctions, we choose some convenient set 
of parameters that clearly display the different features. 
For the tunneling overlap integral we choose the values 
t = 0.2 and t = 0.5 (that correspond, in the notation 
of Ref. to a ~ 0.15 or Z = 2.4 and to a = 0.64 
or Z = 0.75, respectively), and when there is a mag- 
netic field we fix its value to h = 0.2 in units of A (by 
A we mean the magnitude of the singlet gap, Aq, or of 
the triplet vector order parameter depending on the case; 
notice we absorbed Bohr's magneton and the gyromag- 
netic factor in the definition of the magnetic field). These 
values are larger than, for instance, those in the most 
typical STM tunneling experiments, except for the ones 
engineered expressly to seek for large values of but 
have the virtue of making evident the different features, 
including the Andreev gap structure (see below). We 
show, except when indicated, curves for the dc response 
in the limit of vanishing temperatures. For the trunca- 
tion procedure we have taken Na — 3 and TVa = 5 for 
the cases of t = 0.2 and t = 0.5, respectively (and verified 
that larger values produce, given the set of parameters 
chosen, identical curves). The discretization used on the 
horizontal axis is better than 5V — 0.025 A/e in all the 
cases. 

We review now the different pairing-symmetry sce- 
narios. Let us start with the case of normal-metaU 
superconductor junctions. We show in Fig. EKa) typical 
curves for an N-S junction (i.e. a point-contact junction 
between a normal metal and a conventional singlet-paring 
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FIG. 2: Zero temperature I-V characteristics of normal- 
superconductor junctions for both spin-singlet and spin- 
triplet paring, (a) N-S junctions for t = 0.2 (lower curves) 
with and without applied magnetic field (dashed and solid 
line respectively, h — 0.2) and t = 0.5 (upper curve, solid line 
only since the effect of field is smaller and not displayed) ; the 
curves are vertically displaced for clarity, (b) N-T junctions 
for t = 0.2 (upper curve) with and without applied magnetic 
field (dashed and solid line respectively, h = 0.2) and t — 0.5 
(lower curve, solid line only). 
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FIG. 3: Zero temperature I-V characteristics of S-T junctions 
with and without magnetic field (dashed and solid lines re- 
spectively, h = 0.2). The lower solid and dashed curves are 
for t — 0.2 (without and with magnetic field) and the upper 
solid curve is for t = 0.5. The dotted lines are (i) the straight 
unitary slope line for the reference N-N characteristics and (ii) 
the S-S characteristics for similar-parameters junctions (up- 
per t = 0.5 and lower t — 0.2). The inset shows finite temper- 
ature characteristics of S-S junctions with different values of 
the left and right gap amplitudes. The dotted line is the ref- 
erence N-N characteristics and the dashed line is the curve for 
A2 = Ai. The solid lines correspond to A2 7^ Ai (A2 = 3 Ai 
for the line closest to the dashed one and A2 = 37/3 Ai for 
the other one). 



superconductor). The diagonal straight line is the N-N 
characteristics given as a reference. The solid lines corre- 
spond to the N-S junction in zero field and the dashed line 
is for one of the junctions (the less transparent one) in the 
presence of a magnetic field. The effect of the magnetic 
field is to produce what would be seen as a Zeeman split- 
ting of the differential conductance peak (i.e. the peak in 
the curve of dl/dV vs. V). Notice the sub-gap shoulder 
on the I-V curve when eV < A (for instance in the zero 
field case); its origin is in the coherent Andreev processes 
that take place at the junction contact. Next we show 
in Fig. [2J(b) a typical curve this time for what we call 
an N-T junction (i.e. a junction between a normal metal 
and an unconventional triplet-pairing superconductor). 
The solid lines correspond to the N-T junction in zero 
field and the dashed line is for the t = 0.2 junction when 
in the presence of a magnetic field that is aligned with 
the vector order parameter A. If one considers a mag- 
netic field that is perpendicular to the order parameter 
(h _L A), one finds it has no effect on the I-V characteris- 
tic that remains identical to the one for the zero field case 
(remark that in the case of the N-S junction the orienta- 
tion of the field was immaterial) . Notice also the absence 
of a sub-gap shoulder on the I-V curve. This absence is 
caused by the odd real-space symmetry of the supercon- 



ductor (p-wave paring) : Andreev processes with opposite 
chiralities interfere destructively and exactly cancel each 
other. As a result, the curves are exactly identical to 
those computed with a semiconducting band model that 
ignores Andreev scattering (to be contrasted with the 
non-trivial results in this respect that will be shown mo- 
mentarily for junctions involving two different-symmetry 
superconductors) . 

Let us now turn to examine the case of junctions in 
which both their sides are superconducting. In Fig. [3] we 
display typical curves for S-S junctions (both the sides 
are conventional spin-singlet superconductors) and S-T 
junctions (one of the sides is a spin-triplet superconduc- 
tor). The straight dotted line is the N-N characteristics 
- taken as a reference, same as before. The remaining 
dotted lines are the TV curves of S-S junctions that show 
all the standard features already well documented in the 
literature^2ii& For the purpose of later comparison, we 
remark here the sizeable currents for voltages eV > 2A 
(the value of the gap is taken to be the same on both 
sides of the junction), and the 'sub-gap' shoulder with 
Andreev steps at eV = 2A/n (with n = 1,2,3, . . .). We 
also remind the reader that this curve is, when orbital 
effects can be ignored, not sensitive to applied magnetic 
fields. The remaining curves (dashed and solid lines) cor- 
respond to S-T junctions with different tunneling matrix 
element strengths and with and without magnetic field 
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(respectively). The solid lines are insensitive to the orien- 
tation of the vector order parameter on the triple-pairing 
side of the junctions, and the current amplitude is found 
to be systematically smaller than in the case of the re- 
spective S-S junctions. Remarkably, the 'sub-gap' struc- 
ture shows only two steps (at voltages given by n = 1,2) 
and the current is zero when eV < A (if the magnitudes 
of the gaps in the spin-singlet and spin-triplet sides of 
the junction are different, then the zero current condi- 
tion is eV < Axripietj where A Triplet is the magnitude of 
the vector order parameter on the spin-triplet side of the 
junction). Concerning the effects of an applied magnetic 
field, the curves remain unchanged if the field is applied 
parallel to the direction of the vector-order-parameter, 
but show instead a Zeeman effect if the field is perpen- 
dicular to it (dashed line). This is in contrast with the 
case of N-T junctions, for which the Zeeman effect is ex- 
pected for fields h || A. 

In the figure inset we display curves for S-S junctions 
at finite temperatures. In order to render the different 
features simultaneously visible, we 'push' the tempera- 
ture to be equal to A (with = 1). Let us denote as 
A the average gap value (2 A = Ai + A2) that we keep 
using to define the unit in which we measure the voltage, 
normalized current, etc. The dashed line is the finite 
temperature TV for the t — 0.2 junction between two 
identical superconductors (A 2 = Ai), whereas the solid 
lines correspond to similar junctions with A2 = 3 Ai 
or A2 = 37/3 Ai and the same transparency (during 
this discussion we assume A2 > Ai). Besides the stan- 
dard quasiparticle tunneling threshold (eV = 2A) one 
clearly sees in both solid curves the step at eV — A 2 
corresponding to the first term of the so-called even se- 
ries (eV = 2A e /m with £ = 1,2 and m = 2,4, .. .). The 
steps corresponding to the odd series (eV = 2A/m with 
to = 3, 5, . . .) are not visible in this plot, but we verified 
that we observe them at low temperatures in junctions 
with A2 > Ai. This current step structure lacks temper- 
ature or magnetic field dependence (in accordance with 
experimental observations). The other prominent feature 
of the solid curves is the rounded cusp at eV ~ A 2 — Ai. 
This is a thermally activated feature that appears when 
the upper gap edges at both sides of the junctions are 
aligned. The rounding of the cusp has similar origin as 
the rounding of the quasiparticle threshold, both are due 
to higher order multiparticle processes. One new obser- 
vation that we made is that the position of the thermal 
cusp is not exactly A2 — Ai, but it is shifted towards 
lower voltages. This is again a result of taking into ac- 
count higher order Andreev processes and is made more 
evident by our choice of parameters; for lower temper- 
atures and less transparent junctions this correction is 
typically very small. We also find that the steps of both 
the even and the odd series that fall to the left of the 
thermal cusp are usually washed away by its tail; this we 
find to be consistent with available experimental results. 
We expect the different detailed features corresponding 
to dissimilar gaps and finite temperatures to be accessible 
to current state-of-the-art experiments, we will comment 
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FIG. 4: Zero temperature I-V characteristics of T-T junc- 
tions. The lower curve corresponds to t — 0.5 and the upper 
one to t — 0.2. The dotted lines is the reference N-N char- 
acteristics. The inset shows finite temperature characteristics 
of T-T junctions with different values of the left and right gap 
amplitudes. We used the same values choices as in the inset 
of Fig. El Namely, the dotted line is the N-N characteristics 
and the dashed line corresponds to A2 = Ai while the solid 
lines are for A2 7^ Ai(A2 = 3 Ai for the line closest to the 
dashed one and A2 = 37/3 Ai for the other one). 



on that in the next section. 

As an aside, let us comment on the effect of local con- 
tact potential terms in the tunneling matrix. They sup- 
press uniformly the dc current amplitude but have no 
other effect on the shape of the I-V characteristics. This 
is true regardless of the pairing symmetry of the super- 
conductors forming the junction, in particular they do 
not cause the appearance of mid-gap states in the case of 
triplet pairing (for neither the N-T or S-T junctions nor 
for the T-T junctions discussed below). 

Finally, the only remaining case to consider is that of 
junctions in which both sides are spin-triplet supercon- 
ductors. Such a case is exemplified in the curves of Fig.0] 
Andreev processes with triplet symmetric paring, for tun- 
neling through a single- mode contact, interfere destruc- 
tively and the current remains zero up to voltages larger 
than eV = 2A, when quasiparticle tunneling becomes al- 
lowed. The lower solid line corresponds to t = 0.5 and the 
upper one to t = 0.2, the inversion of the order is due to 
the fact that, for this range of parameters, the current at 
fixed voltage grows more slowly than in the case of nor- 
mal junctions (N-N). Similarly as in the S-S case, both 
sides of the junction react identically to applied mag- 
netic fields and no net effects are therefore visible in the 
current- voltage characteristics. 

The curves in the figure inset are the same as in the 
inset of the previous figure but for the case of spin-triple 
pairing symmetry (on both sides of the junction). Both 
the quasiparticle tunneling threshold and the thermal 
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cusp remain sharp since higher order processes interfere 
destructively and no rounding takes place. Another con- 
sequence of this is that the position of the thermal cusp 
is exactly eV = A 2 — Ai and no shift is observed. The 
sub-gap part of the curves is smooth and its height and 
shape are governed by the thermal excitations, the step 
structure of the even and odd series is absent. Also in 
this case magnetic fields have no direct effects. 



V. EXPERIMENTAL CONSEQUENCES 

A first set of applications concerns atomic contacts. 
For the case of normal or singlet superconducting leads 
with identical gaps, and zero temperature and magnetic 
field, our results are in full agreement with the previous 
studies of such systems, 16 For the case of two different 
gaps shown in the inset of Fig. our theory correctly re- 
produces the different steps as discussed in the previous 
section. Another prominent feature of such curves is the 
thermally activated rounded cusp at eV « A2 — Ai. The 
rounding is due to higher order multiparticle processes; 
and a new prediction is that the position of the thermal 
cusp is not exactly A2 — Ai, but it is shifted towards 
lower voltages. Such a shift from the naive "density of 
states" answer could in principle be checked directly in 
atomic contacts. We also find that the steps of both the 
even and the odd series that fall to the left of the thermal 
cusp are usually washed away by its tail; this we find to 
be consistent with available experimental resultsi Such 
features corresponding to dissimilar gaps and finite tem- 
peratures should be accessible to current state-of-the-art 
experiments. For instance, the experiments of Ref. [l^ 
could be attempted using a Pb STM tip as before but to 
prove into a Mn-doped Pb sample. Mn will act as a mag- 
netic impurity and decrease the value of the gap, such a 
setup would correspond to the situation of similar but 
not identical gap parameters with a doping controllable 
difference; it would be a way to try to observe the 'split- 
ting' of the even series in single-point contacts. Other 
setups could be envisaged based also on STM techniques 
or on pressed crossed wires. 

The main application of our results, however, con- 
cerns the use of tunneling with triplet superconductors. 43 
In that case the most direct experimental realization is 
organic superconductors.^ The experiments show that 
for magnetic fields along the direction of the conducting 
chains (a crystalline-axis) the upper critical field is para- 
magnetically limited. If such systems are indeed triplet 
superconductors, this would correspond, following our 
notations, to a vector order parameter aligned with the 
field (h || A)iSi With this geometry a Zeeman splitting of 
the differential conductance peak, similar to that in con- 
ventional superconductors, should be observed in a tun- 
neling experiment. As the field is rotated the splitting 
would be suppressed and for a magnetic field oriented 
parallel to the b' crystalline- axis there should be no Zee- 
man effect (accompanied by the possibility of applying 



large fields that are not paramagnetically limited). The 
disappearance of splitting even as the field is being in- 
creased would constitute a clear signature of spin-triplet 
superconductivity. The main difficulties of such an ex- 
periment would be the set up of point-contacts and the 
resolution required to observe the Zeeman effect. On the 
first point one possibility would be to use STM setups 
with 'thin tips'. On the second point, since the criti- 
cal temperature of these organic salts is relatively low, 
the experiment could be done with moderate fields that 
would produce splittings that are a substantial fraction 
of the superconducting gap. The linearity of the mag- 
netic field dependence of these splittings, a signature of 
the Zeeman effect, could be accurately established using 
Fourier analysis techniques^ 

Similarly as in the case of N-T junctions, we can en- 
visage using the Zeeman response of S-T junctions as 
a direct probe for spin-triplet order. If, for instance, 
a magnetic field is applied along the b' crystalline-axis 
of (TMTSF) 2 PF 6 , we predict a Zeeman splitting of the 
main differential conductance peak. This would also con- 
stitute a clear sign of unconventional superconductivity 
since such an effect does not take place for standard BCS 
superconductors. The b' direction is the one on which 
the upper critical field is not paramagnetically limited, so 
relatively large fields could be applied in order to obtain 
a clear signal (as the field alignment changes the split- 
ting disappears) . To afford large fields one would need to 
use in the 'conventional' side of the junction a compound 
with a relatively high critical temperature (as compared 
with that of (TMTSF) 2 PF 6 ). In this respect one has the 
bonus that, since the required setup should be a point- 
contact, superconductivity might survive at the contact- 
neck region up to fields much in excess of the bulk value 
of H C 2 (rather approaching Hp CS for that material) 
Another advantage in the two-superconductor setup is 
that the levels of noise are usually smaller)- 3 allowing 
a better definition of the differential conductance signal 
from where the Zeeman splitting is going to be read off. 

Similar considerations could also be made for those 
layered compounds that are believed to be triplet 
superconductorsiS*^ In that case the critical magnetic 
field is not paramagnetically limited when the applied 
field is oriented parallel to the superconducting planes. 
Among these compounds Sr2Ru04 is the best studied 
one so far, but only few tunneling experiments were 
performed r^S^LSiS and none so far with high resolu- 
tion and in the presence of an applied external magnetic 
field (see though Refs. |57ll58ll60| ). One of the conspic- 
uous features observed in some of these experiments is 
the presence of a 'zero bias anomaly' (ZBA) in the dif- 
ferential conductance. Its explanation is still a matter of 
debate, but seems to require extended contact interfaces 
and momentum dependent order parameters (to include 
the effect of 'zero energy states' at the interfaces, exten- 
sions to our scheme would be required, possibly incorpo- 
rating certain aspects of those calculations already done 
for planar junctiona^i*!^ 3 .). Our general findings about 
the effect of magnetic fields should, however, apply, since 
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they refer to effects to be measured at voltages of the 
order of the superconducting gap. It is intriguing to no- 
tice that in the 'point-contact' experiment of Ref. |Hfij 
two types of spectra are measured: with and without 
ZBA in the differential resistance. One might speculate 
that what changes between the different samples should 
be no other thing than the effective size, potential bar- 
rier and geometry of the contact (experiments in other 
compounds indicate that that might be enough to give 
rise to zero voltage features; cf. Ref. H3). In particular, 
if that is the case, at least those dV/dl curves with no 
ZBA should, according to our calculations, show no Zee- 
man splitting when a held is applied parallel to the Ru-0 
planes; in contrast to what is expected for BCS supercon- 
ductors. Two different groups reported that further point 
contact and STM tunneling experiments on Sr2Ru04 in 
a magnetic held are underway^S^ 

VI. SUMMARY AND CLOSING REMARKS 

Summarizing, we have shown how the full I-V char- 
acteristics for point-contact junctions can be accurately 
studied using a local action approach in the context of the 
Keldysh formalism. Our formalism allows one to treat 
both normal and superconducting (singlet and triplet) 
leads, and to take into account effects of finite magnetic 
field and temperature. In particular we have shown that 
the point-contact tunneling involving unconventional su- 
perconductors with spin-triplet pairing displays interest- 
ing characteristic features. Unlike the case of conven- 
tional superconductors, these show quite different char- 
acteristics whether the junctions are planas^ or point- 
contact like42i The Zeeman response to an external mag- 
netic field is such that it allows for the identification of 
triplet phases and might be relevant for future experi- 
ments. The prediction of a truncated sub-gap structure 
in point-contact S-T junctions is also very interesting, 
but experiments to test this are much harder to carry 
out. That kind of detailed experiments are, however, 
possible for conventional superconductors and we believe 
they could be exploited to look for some as yet poorly 
tested predictions of the theory. For instance, the exper- 



iments of Ref. [3] could be attempted using a Pb STM 
tip as before but doping the sample with Mn; it would be 
a way to try to observe the 'splitting' of the even series 
in single-point contacts. 

Besides the different additional effects on the tunnel- 
ing characteristic that we discussed here (the effects of 
fields, temperature and contact potentials), there are 
others that can also be easily taken into account like, 
for instance, spin-flip tunneling processes or temperature 
gradients. These effects will be relevant in the study of 
junctions involving ferromagnets (of possible relevance 
in the context of spintronics; cf. Ref. or in preci- 

sion studies related to the renewed interest in the use of 
microjunctions for on-chip thermometry and eventually 
cryogenics*^ A remaining challenge, however, is how to 
extend our formalism seeking to include the physics re- 
sponsible for producing ZBAs in order to further our un- 
derstanding of tunneling experiments in layered uncon- 
ventional superconductors and planar junctions. Such 
extensions shall seek to describe not only the point- 
contact limit, but also the planar interface case. This is 
one of the ingredients necessary for a realistic description 
of layered materials like the ruthenates. In that respect 
the interplay with the advances already made using semi- 
classical methods will constitute not only a check of the 
approximations made in the latter but also a way of find- 
ing efficient computational schemes for more complicated 
scenarios that might incorporate, for instance, the two- 
band nature of certain compounds. On the other hand, 
in the context of the quasi-one-dimensional organic con- 
ductors, such developments should help to go beyond the 
one-dimensional approximation. 

Acknowledgments 

We would like to thank 0. Fischer, M. Eskildsen, 
M. Kugler and G. Levy for discussions about tunneling 
and STM. We also thank Y. Maeno and M. Sigrist for 
discussions about tunneling into triplet superconductors 
and ruthenates in particular. This work was done under 
the auspices of the Swiss National Science Foundation 
and its MaNEP program. 



1 E. L. Wolf, Principles of Electron Tunneling Spectroscopy, 
vol. 71 of International Series of Monographs on Physics 
(Oxford University Press, Clarendon Press, Oxford, 1989), 
2nd ed. 

2 C. F. Quate, Phys. Today 39, 26 (1986). 

3 G. Binnig and H. Rohrer, Rev. Mod. Phys. 59, 615 (1987). 

4 G. Binnig and H. Rohrer, Rev. Mod. Phys. 71, S324 
(1999). 

5 J. Tersoff and D. R. Hamann, Phys. Rev. Lett. 50, 1998 
(1983). 

6 J. Tersoff and D. R. Hamann, Phys. Rev. B 31, 805 (1985). 

7 H. Suderow, E. Bascones, W. Belzig, F. Guinea, and 
S. Vieira, Europhys. Lett. 50, 749 (2000). 

8 H. Suderow, E. Bascones, A. Izquierdo, F. Guinea, and 



S. Vieira, Phys. Rev. B 65, 100519(R) (2002). 
9 J. C. Cuevas, A. L. Yeyati, and A. Martm-Rodero, Phys. 
Rev. Lett. 80, 1066 (1998). 

10 J. Nicol, S. Shapiro, and P. H. Smith, Phys. Rev. Lett. 5, 
461 (1960). 

11 T. M. Klapwijk, G. E. Blonder, and M. Tinkham, Physica 
B+C 109-110, 1657 (1982). 

12 G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 
Rev. B 25, 4515 (1982). 

13 M. Octavio, M. Tinkham, G. E. Blonder, and T. M. Klap- 
wijk, Phys. Rev. B 27, 6739 (1983). 

14 M. H. Cohen, L. M. Falicov, and J. C. Phillips, Phys. Rev. 
Lett. 8, 316 (1962). 

1 J. W. Wilkins, Tunneling Phenomena in Solids (Plenum 



10 



Press, New York, 1969), chapter 24, p. 333, Multiparticle 
Tunneling, 1st ed. 

16 J. C. Cuevas, A. Martm-Rodero, and A. L. Yeyati, Phys. 
Rev. B 54, 7366 (1996). 

17 E. Scheer, P. Joyez, D. Esteve, C. Urbina, and M. H. De- 
voret, Phys. Rev. Lett. 78, 3535 (1997). 

18 E. Scheer, N. Agrai't, J. C. Cuevas, A. L. Yeyati, B. Lu- 
doph, A. Martm-Rodero, G. Rubio-Bollinger, J. M. van 
Ruitenbeek, and C. Urbina, Nature 394, 154 (1998). 

19 B. Ludoph, N. van der Post, E. N. Bratus', E. V. Bezuglyi, 
V. S. Shumeiko, G. Wendin, and J. M. van Ruitenbeek, 
Phys. Rev. B 61, 8561 (2000). 

20 E. Scheer, W. Belzig, Y. Naveh, M. H. Devoret, D. Esteve, 
and C. Urbina, Phys. Rev. Lett. 86, 284 (2001). 

21 G. Rubio-Bollinger, C. de las Heras, E. Bascones, 
N. Agrait, F. Guinea, and S. Vieira, Phys. Rev. B 67, 
121407(R) (2003). 

22 M. Hafner, P. Konrad, F. Pauly, J. C. Cuevas, and 
E. Scheer, Conduction channels of one-atom zinc contacts, 
arXiv:cond-mat /0407207. 

23 A. Millis, D. Rainer, and J. A. Sauls, Phys. Rev. B 38, 
4504 (1988). 

24 J. C. Cuevas and M. Fogelstrom, Phys. Rev. B 64, 104502 
(2001). 

25 J. Kopu, M. Eschrig, J. C. Cuevas, and M. Fogelstrom, 
Phys. Rev. B 69, 094501 (2004). 

26 N. Kopnin, Theory of Nonequilibrium Superconductivity, 
vol. 110 of International Series of Monographs on Physics 
(Oxford Science Publications, Clarendon Press, Oxford, 
2001), 1st ed. 

27 P. W. Anderson and P. Morel, Phys. Rev. 123, 1911 
(1961). 

28 R. Balian and N. R. Werthamer, Phys. Rev. 131, 1553 
(1963). 

29 D. D. Osheroff, R. C. Richardson, and D. M. Lee, Phys. 
Rev. Lett. 28, 885 (1972). 

30 D. D. Osheroff, W. J. Gully, R. C. Richardson, and D. M. 
Lee, Phys. Rev. Lett. 29, 920 (1972). 

31 P. W. Anderson and W. F. Brinkman, in The Helium Liq- 
uid, edited by J. G. M. Armitage and I. E. Farquhar (Aca- 
demic Press, New York, 1975), Proceedings of the 15 th 
Scottish Universities Summer School in Physics, pp. 314- 
416. 

32 Chem. Rev. (2004), special issue on Molecular Conductors. 

33 T. fshiguro, in High Magnetic Fields: Applications in 
Condensed Matter Physics and Spectroscopy, edited by 

C. Berthier, L. P. Levy, and G. Martinez (Springer- Verlag, 
Heidelberg, 2002), vol. 595 of Lecture Notes in Physics, pp. 
301-313. 

34 T. M. Rice and M. Sigrist, J. Phys.: Condens. Matter 7, 
L643 (1995). 

35 G. Baskaran, Physica B 224, 490 (1996), proceedings of 
the International Conference on Strongly Correlated Elec- 
tron Systems. 

36 R. Joynt and L. Taillefer, Rev. Mod. Phys. 74, 235 (2002). 

37 J. de Phys. IV 114 (2004), ISCOM 2003 - The Fifth Inter- 
national Symposium on Crystalline Organic Metals, Super- 
conductors and Ferromagnets. Port-Bourgenay, September 
21 s '-26'\ 2003. 

38 J. I. Oh and M. J. Naughton, Phys. Rev. Lett. 92, 067001 
(2004). 

39 N. Joo, P. Auban-Senzier, C. R. Pasquier, P. Monod, 

D. Jrome, and K. Bechgaard, Eur. Phys. J. B 40, 43 
(2004). 

40 T. M. Rice and M. Sigrist, Phys. Today 54, 42 (2001). 

41 A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 



(2003). 

H. I. Ha, J. I. Oh, J. Moser, and M. J. Naughton, Synth. 
Metal 137, 1215 (2003). 

C. J. Bolech and T. Giamarchi, Phys. Rev. Lett. 92, 
127001 (2004). 

I. Affleck, J. S. Caux, and A. M. Zagoskin, Phys. Rev. B 
62, 1433 (2000). 

T. Giamarchi, Quantum Physics in One Dimension, vol. 
121 of International Series of Monographs on Physics (Ox- 
ford Science Publications, Clarendon Press, Oxford, 2004), 
1st ed. 

G. D. Mahan, Many-Particle Physics, Physics of Solids and 
Liquids (Kluwer Academic/Plenum Publishers, New York, 
2000), 3rd ed. 

L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965). 

Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 

(1999). 

C. Grimaldi and P. Fulde, Phys. Rev. B 56, 2751 (1997). 
A. G. Lebed, K. Machida, and M. Ozaki, Phys. Rev. B 62, 
R795 (2000). 

G. B. Hertel and T. P. Orlando, Phys. Rev. B 32, 166 
(1985). 

A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 
M. Eskildsen, private communication. 

V. Mineev, in High Magnetic Fields: Applications in 
Condensed Matter Physics and Spectroscopy, edited by 
C. Berthier, L. P. Levy, and G. Martinez (Springer- Verlag, 
Heidelberg, 2002), vol. 595 of Lecture Notes in Physics, pp. 
287-300. 

R. Jin, Y. Liu, Z. Q. Mao, and Y. Maeno, Europhys. Lett. 
51, 341 (2000). 

F. Laube, G. Goll, H. v. Lohneysen, M. Fogelstrom, and 

F. Lichtenberg, Phys. Rev. Lett. 84, 1595 (2000). 

Z. Q. Mao, K. D. Nelson, R. Jin, Y. Liu, and Y. Maeno, 
Phys. Rev. Lett. 87, 037003 (2001). 

M. D. Upward, L. P. Kouwenhoven, A. F. Morpurgo, 
N. Kikugawa, Z. Q. Mao, and Y. Maeno, Phys. Rev. B 
65, 220512(R) (2002). 
M. Kugler, private communication. 

A. Sumiyama, T. Endo, Y. Oda, Y. Yoshida, A. Mukai, 
A. Ono, and Y. Onuki, Physica C 367, 129 (2002). 
M. Yamashiro, Y. Tanaka, and S. Kashiwaya, Phys. Rev. 
B 56, 7847 (1997). 

C. Honerkamp and M. Sigrist, J. Low Temp. Phys 111, 
895 (1998). 

Y. Asano, Y. Tanaka, M. Sigrist, and S. Kashiwaya, Phys. 
Rev. B 67, 184505 (2003). 

G. Sheet, S. Mukhopadhyay, and P. Raychaudhuri, Phys. 
Rev. B 69, 134507 (2004). 

H. Kambara, K. Yokota, T. Matsui, I. Ueda, T. Shishido, 
M. Wada, N. Kikugawa, Y. Maeno, and H. Fukuyama, 
Physica C 388-389, 503 (2003), proceedings of the 23 rd 
International Conference on Low Temperature Physics 
(LT23), Hiroshima, Japan, 20-27 August 2002. 

K. Sengupta, I. Zutic, H.-J. Kwon, V. M. Yakovenko, and 

S. Das Sarma, Phys. Rev. B 63, 144531 (2001). 

R. Lopez and D. Sanchez, Phys. Rev. Lett. 90, 116602 

(2003). 

J. Pekola, R. Schoelkopf, and J. Ullom, Phys. Today 57, 
41 (2004). 

Notice that the object aipLa transforms as the complex 
conjugate representation of the fundamental representa- 
tion of SU(2). In other words, it has the same transforma- 
tion properties as V4, a - 



